SPONTANEOUS VS. PIEZOELECTRIC POLARIZATION IN III-V NITRIDES: 
CONCEPTUAL ASPECTS AND PRACTICAL CONSEQUENCES 
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Macroscopic polarization plays a major role in determining the optical and electrical properties of 
nitride nanostructures via polarization-induced built-in electrostatic fields. While currently fashion- 
able, this field of endeavour is still by far in its early infancy. Here we contribute some clarifications 
on the conceptual issues involved in determining built-in fields in III-V nitride nanostructures, sort- 
ing out in particular the roles of spontaneous and piezoelectric polarization. 

Subject classification: 73.90+f,77.65-j,77.90+k,S7.14 



o 

CZ5 



-a 
c 

o 
o 



> 

00 

o 

00 

o 

On 
On 
-i— > 



i 

C 

o 
o 



I. INTRODUCTION 

Macroscopic polarization in low-symmetry crystals is 
well known in the ferroelectrics community. Outside 
that community, it was so far held for a curiosity de- 
void of practical, let alone technological, relevance. In 
particular, it was known that a large spontaneous po- 
larization (henceforth SP), at most an order of magni- 
tude smaller than the giant values of ferroelectric per- 
ovskites, was present in some wurtzite semiconductors, 
but due to their minor technological relevance compared 
to zincblende III-V arsenides or phosphides, this issue 
was largely overlooked outside a limited circle of special- 
ists. 

Such state of affairs has changed after first-principles 
calculations || provided evidence of the considerable SP 
and larger-then-usual piezoelectric coupling constants of 
wurtzite III- V nitrides. Because of the huge importance 
of the latter as materials for optoelectronic devices in the 
blue-UV range H , the issue of spontaneous polarization 
fields was brought to a larger community. 

The obvious relevance of polarization in nitrides is due 
to the fact that it induces large (~MV/cm) built-in elec- 
trostatic fields in layered nitride nanostructures, as first 
shown in Ref . ^| . These fields affect dramatically the op- 
tical and electrical properties of those structures [5-14]. 
Quite a flurry of activity has recently revealed a host of 
new and unexpected effects: to name but some, Stark- 
like red shifts in recombination and absorption energies 
for increasing quantum well widths p|~p| , concurrent sup- 
pression of oscillator strength pM pfl, ensuing anomalies 
in recombination dynamics |p|]lC|-|l^j, and their interplay 
with charge injection [||^[l3| driving a recovery of oscilla- 
tor strength and a blue back-shift of transition energies, 
"self-doping" effects at HEMT heterointerfaces |l4|] , have 
been demonstrated theoretically and experimentally. 

In view of these exciting developments, it is likely that 
most of the potentially observable polarization-related ef- 
fects in nitrides are still to be dreamt about. It is safe 
to state, in any event, that conceptual clarity on the 
determination of the electric fields produced by polar- 



ization, and in particular by SP, is a must if we are to 
understand and exploit these effects in practical appli- 
cations. To mention just a pair of issues, the built-in 
fields depend on boundary conditions and on device de- 
sign (number of quantum wells, thickness, composition, 
contacts, etc.); also, the field values and sign patterns 
are not those expected if only piezoelectricity is consid- 
ered. This work therefore aims at clarifying some of these 
subtle questions, focusing on three inter- related systems, 
i.e. large samples, multi-quantum- wells (MQWs), and 
isolated quantum- wells (QWs). 



II. SPONTANEOUS VS. PIEZOELECTRIC 

Spontaneous polarization and piezoelectric constants 
can be reliably determined [15| via simple electronic 
structure calculations based on the Berry geometric 
quantum phase concept. We start by discussing our re- 
cent result thereon [0], which are summarized in Figs. 
|l| and ||. The values calculated ab initio are for III-N 
binaries (A1N, GaN, InN) only; since direct calculations 
for the alloys are not yet available, an estimate of the 



SP P( sp > in Al x In J/ Ga (1 _ x _j /) N alloys was obtained by 
Vegard interpolation: 
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similar relations holding for the piezoelectric constants. 
Thus, the triangle in Fig. |l| borders the (theoretical) SP 
values achievable in a general III-N alloy as a function of 
in-plane lattice constant. The triangles in Fig. || border 
the accessible values of the two piezoconstants relevant 
to a-plane epitaxial strains. The large range of accessible 
polarizations and lattice constants is due to the concur- 
rent large lattice mismatch between GaN and InN on one 
side, and the very large SP in A1N on the other. 

Three main points are evident from the data reported. 
First, the piezoelectric constants are an order of mag- 
nitude larger than in other III-V's, and reversed in 
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sign. Second, for basal-plane strains typical of epitax- 
ial wurtzite nitride multilayers (2-5%), SP is comparable 
to, or larger than piezoelectric polarization. Therefore, 
SP is all but negligible in any typical III-V nitride nanos- 
tructures. Third, GaN has nearly the same SP as InN, 
but a large lattice mismatch to it; on the other hand, 
A1N has a SP about three times larger than GaN, but a 
much smaller mismatch. This implies that InGaN/GaN 
structures will tend to be mostly influenced by piezoelec- 
tricity, whereas in AlGaN/GaN systems, SP effects will 
be dominant. Mixed regimes can also occur, which we 
will come back to in Sec. M. 



III. MASSIVE SAMPLES 

By elementary electrostatics, polarization and electro- 
static field are related by 

E=(D-P)/ £ (2) 

where D is the displacement field, P the macroscopic po- 
larization, and e the static dielectric constant. The latter 
have been calculated to be 10.31, 10.28, and 14.61 
for A1N, GaN, and InN respectively; the value for generic 
alloys are estimated by Vegard interpolation. The polar- 
ization in Eq. || is known as transverse, and it is the sum 
of the spontaneous and piezoelectric polarization contri- 
butions: P = P( s p) + P(p z ). The displacement field is 
determined by the free-charge distribution in the mate- 
rial: 

V-D = e(p-ra), (3) 

with e the electron charge, and p and n the hole and elec- 
tron densities, respectively. The evaluation of the electric 
field E inside the semiconductor requires in general a self- 
consistent solution of Eq. and ||| as outlined in Refs. || 
and |^. However, many cases can be discussed without 
explicit calculations, an useful exercise especially in the 
low free-carrier-density limit. 

An interesting such case is the evaluation of polariza- 
tion effects in large, massive samples. In effect, misun- 
derstanding the role of the displacement field in massive 
samples leads to predicting wrong electrostatic fields in 
a MQWs system. At the outset, let it be pointed out 
that a strong polarization field in a III-V nitride massive 
sample does not imply the existence of a macroscopic 
electrostatic field over the whole sample. Qualitatively, 
since such an electric field would produce (at non-zero 
temperature) a persistent current due to thermally gen- 
erated intrinsic carriers, and because this current cannot 
be sustained indefinitely, the electric field inside a mas- 
sive sample must be zero. One deduces thereforth that 
the displacement field in the sample is 

D = P (s P )_ (4) 



For a neutral (albeit polarized) sample, displacement 
field conservation across the surface would imply an elec- 
tric field E = P( s p) in the vacuum off the sample surface. 
Of course this is unrealistic, as outside the sample P, 
E, and D must all vanish. The SP discontinuity across 
the surface must be counterbalanced by an equal change 
in the displacement field. Thus, near the surface of our 
sample, there exists a region where the free-carriers dis- 
tribution satisfies 

J (e(p-n)-V-P (sp) )dz = 0, (5) 

with £ the thickness of the surface region, which may 
range easily in the hundreds of A. In this region 
the free carriers will constitute an approximately two- 
dimensional electron gas (2DEG) whose areal density 
equals the SP. The electric field in this region depends 
strongly on material parameters and on surface states. It 
can be concluded however that at distances in excess of £ 
from the surface (or, in the case of devices, the interface 
with contacts, buffers, caps, etc.), the electric field will 
be zero, and the polarization will equal p( s P), as will the 
displacement field. This is the starting point to discuss 
MQWs. 

IV. MULTI-QUANTUM- WELLS 

As pointed out above (especially in a semi-insulating 
sample) it can be safely assumed that the displacement 
field in a large sample is uniform all over the crystal, and 
equal to p( sp ), and E = 0. If a sequence of QWs, of fixed 
composition, is now inserted in the (otherwise homoge- 
neous) sample, the electric field inside the generic j-th 
layer (either QW or barrier) is given by: 

E 3 = (D-P 3 )/ £j (6) 

with Pj and Ej the total polarization and dielectric con- 
stant in layer j. 

If the displacement field is constant, i.e. if Eq. || holds, 
the electric field in the j-th layer is 

% = (pW-P J .)/ £i) (7) 

where p( sp ) is the SP in the massive sample material, 
which now functions locally as barrier material. The 
field is then zero in the barriers, which are made up of 
unstrained material with spontaneous polarization p( s P), 
while each QW in the series is subject to the same electric 
field Ej. A potential drop given by 

= |Ej| • lj (8) 

will thus occur in the generic j-th well of thickness lj. 
Now Eq. [I] holds only subject to the following four con- 
straints: 

i) the Debye-Hiickel screening length is at least larger 
than, say, the MQW region thickness; 
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ii) the QWs are far enough from any interface with the 
outer world that the surface effects mentioned in the pre- 
vious Section are negligible; 

Hi) none of the AVi, nor their sum, must exceed the band 
gap energy of the well material: 



^ ZfeEfe < E g 



(9) 



can be calculated by simply solving a Schrodinger equa- 
tion for the compositional potential of the superlattice 
with the built-in fields as given by Eq. [l2| Of course, 
this is emphatically not the case in general, because of 
e.g. doping, excitation, boundary conditions.etc. A self- 
consistent solution as outlined e.g. in Refs. g and ^ will 
be needed. 



iv) the well-barrier interfaces are free of interface states. 

These requirements are typically met in nitride 
MQWs, except for condition (Hi), which breaks down if 
the QWs form a sufficiently thick superlattice. Thereby, 
a self-consistent calculation explicitly accounting for free 
carriers must be performed. To approximately enforce 
condition (Hi) avoiding cumbersome calculations, it is 
convenient to approximate Eq. ^| applying periodic 
boundary conditions, 



^ fc E fc = 0, 



(10) 



where the sum runs over all layers in the MQW, in par- 
ticular including the barrier layers. Thereby the displace- 
ment field is determined subject to the condition of zero 
average electric field in the MQW region. The maximal 
error in the predicted fields due to using Eq. |o| instead 
of Eq. |^ is of the order of E gap /d, with d the overall 
thickness of the superlattice; hence this approximation 
becomes exact in the limit of very thick superlattices. 
The error may be of course smaller depending on the 
specific conditions in the device (doping, contacts, etc.). 

Substituting Eq. ^ in Eq. [t(J the displacement field is 
then obtained: 



D 



Ysk l k/ £ k 



(11) 



Plugging this back into Eq. ^ gives a simple expression 
for the field in the generic well or barrier of the MQW, 
namely 



E; 



Efc IkPk/fk - Pj Efc l k/s k 

£ 3 Lfe l k/£k 



(12) 



with sums running on all layers (including the j-th). This 
is a general expression for any thickness of wells and bar- 
riers in a generic superlattice or MQW, which reduces to 
the formulas in Ref. [M] assuming a single type of well 
and barrier. 

It is to be noted again that in this case the fields are 
non-zero both in the wells and the barriers. Typically, 
the sign of the barrier field will be opposite to that of 
the well, leading to an effectively triangular confinement 
potential, which favors transfer of confined electrons into 
the barrier layer on one side of the well, which may lead 
to anomalous effects in excitonic transitions |i"H ]. 

In practice, if the conditions discussed above are sat- 
isfied for a specific system, the properties of the latter 



V. ISOLATED QUANTUM WELLS 

Consider now a single QW of material A; let the latter 
have dielectric constant sa and spontaneous polarization 
P^ p ^ , and the QW made thereof be strained so as to be 
subject to a piezoelectric polarization p( pz ). The total 
polarization in the QW is then P^ p) + P (pz) . The well 
is embedded in an extended, insulating, and unstrained 
material with spontaneous polarization p( sp ). According 
to Eq. §, the field in the QW is 



E = [P (sp) 
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]/ £A -P (pZ) /£A 



(13) 



At variance with conventional III-V's (having no SP), an 
additional term appears, due solely to the differ- 

ence in the SP between the QW active layer and the bar- 
rier material p7| . This additional term has far-reaching 
consequences, and can be used to tune the value of the 
electric field inside the QW. In Fig. ||we depict the value 
of the electric field E and of its purely piezoelectric com- 
ponent E( pz ) vs in-plane strain for an AlInGaN QW em- 
bedded in GaN. Without SP, the values of the electric 
field inside the QW would fall into the hatched region 
delimited by dashed lines. The inclusion of SP gives 
rise to a different, wider accessible region, especially for 
small strains. The region in question is white and delim- 
ited by a continuous line. On the negative strain side of 
Fig. ||, this region barely touches the upper side of the 
piezoelectric region, and shrinks becoming a line for large 
strains. The various region boundaries are curved due to 
the quadratic dependence of the piezoelectric component 
on alloy composition (through [piezoconstants] x [strain] 
terms). 

Several features are worth a mention. First, the two 
regions do not overlap, i.e. for any composition the total 
field differs from its pure piezoelectric component (this of 
course is due to the fact that E^ sp ^ never vanishes). This 
difference is rather dramatic for positive strains, i.e. in 
Al-rich alloy wells. A noticeable exception to this general 
behavior is that of GaN/InGaN QWs. These correspond 
to the region where the total and piezo field regions in 
Fig. H (almost) touch each other: in that case indeed, as 
InN has nearly the same SP as GaN, the SP difference 
between barrier and active layer, and hence E^ sp \ is neg- 
ligible with respect to the piezoelectric term. That is why 
early investigations on GaN/InGaN MQWs [^] yielded 
results in good agreement with theory even though SP 
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was neglected. Second, without SP the electric field is 
zero at zero strain (the typical situation in zincblende 
III-Vs). This is no more the case in III-V nitrides, where 
SP is an additional degree of freedom potentially produc- 
ing fields up to 4 MV/cm at zero strain. This requires 
growing AllnGaN with appropriate compositions, which 
appears to be difficult because of thermodynamic solu- 
bility constraints; however, fields of several hundreds of 
KV/cm can be achieved already at minute Al and In 
concentrations, which may perhaps become accessible in 
the future. Third, last but not least, SP provides a han- 
dle for reducing or even, in principle, zeroing the field in 
strained QWs Specifically, this occurs on the zero 
field line in Fig. [|. This can be a breakthrough for the 
many applications needing a good confinement, but at 
the same time no built-in field in the active region. The 
same considerations on growth hold as in the previous 
point. 



VI. SUMMARY AND ACKNOWLEDGEMENTS 

In conclusion, SP can be considered as a degree of free- 
dom to tune the value of the polarization-induced electric 
fields in QW systems. It was shown, for a GaN/ AllnGaN 
QW, that fields up to 4MV/cm can be obtained in ab- 
sence of strain, and that, conversely, a vanishing field can 
also be obtained, despite lattice mismatch, for judicious 
choices of composition and strain. 
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FIG. 1. Spontaneous polarization in AlaJn^Gai-x-j/N alloys according to a Vegard-like rule. 
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FIG. 2. Piezoelectric coupling constants e33 (squares) and e3i (dots) after Ref. |2| Triangles refer to AlInGaN alloys according 
to a Vegard-like rule. 
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FIG. 3. Electric field in an isolated ALTn^Gai-^-yN QW with GaN barriers as a function of the well in-plane strain for 
all possible x and y. The white region delimited by the solid line encloses the accessible values of the total electric field, the 
hatched region limited by dashed line encloses those of the purely piezoelectric polarization component. The two lines are 
barely in contact at the upper edge of the hatched region: this contact zone, indicated by the arrow, corresponds to the InGaN 
alloy case (see text) . 
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